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\ (Icm  ri|)|:on  of  thr  rtr-lt-vrl  blnoT-prcssnre  rcspjiue  (o  + GZ 
'lri',s  in  rrliixrii  hiimsins  wax  nhlnlrici'  «v  an  empirical,  ensemble- 
awrai’c.  (Mn-hlimil-prr.ssurr  transfe.  function  based  on  three 
MibieiTs  responses  to  a total  of  21  simulated  tlerial  combat 
nianeinerini!  runs  on  (be  t 'SAKS AM  human  centrifuge.  Three 
ilifl irent  analsfic  transfer  functions  ssere  fitted  to  the  empirical 
fom  lion  for  frequencies  front  5 to  2i  0 mil/,  and  predictive  per- 
formance of  the  empirical  and  three  analytic  functions  was  ex- 
amined. I he  double-zero,  dr  uhle-pol  mathematical  model  most 
tinsels  til  the  empirical  transfer  fun  lion  and  displayed  reason- 
able prcdlcli'e  ability. 


A ( ISI-  DESCRIPTION  nf  the  human  cyc-lcvc! 
- V arterial  blood-press  ire  response  to  +Gz  stress  is 
lack ini».  !i  is  well  estii  'lislicd  (2)  that  a hydrostatic 
ei'lumn  efleet  nets  to  educe  blood  pressure  at  eye 
Iced  below  that  at  heart  level:  the  amount  of  this  rc- 
c 1 tic  t ion  would  be  equal  to  cl  * ;i  * I , where  d is  the  density 
of  oloml,  ;•  is  the  ambient  acccleratnry  field,  and  It  is 
die  vertical  heart-to-cye  distance  (approximately  30 
cm),  it  the  cardiovascular  system  were  to  respond  as  a 
passive,  rigid-walled  system,  in  such  a system,  the 
change  in  eye-level  blood  pressure  during  G stress  could 
be  described  simply  as  -3.0  X 104  dyncs/cm'VG,  or 
about  -22  mm  Ilg/G.  .because  the  arterial  system  is  n 
distensible  hydraulic  column,  tis  oointcul  out  by  Lawton, 
cf  al.  (')),  and  the  cardiovascular  (baroceptor)  reflexes 
serve  to  counter  G-stress-indu  cd  reductions  in  cyc- 
level  arterial  pressure  (3,6,1  I),  he  description  of  blood 
pressure  tis  a function  of  G load  and  time  is  consider- 
ably complicated,  Mathematical  models  of  the  relation 
between  O load  and  blood  pressure  have  been  based  on 
animal  experiments  (5,15)  anti  on  deconvolution  of 
scarce,  published  human  blood-pressure  response 
curves  (12),  hut  adequate  data  upon  which  to  base  a 

The  research  reported  ir  this  paper  wns  conducted  by  per- 
sonnel of  the  C i v Technology  and  biometrics  Divisions,  USAP 
School  of  Aeros  ice  Medicine 

J,  J.  I reermu  is  n US/  T;-ASF,R  bummer  Faculty  Research 
Program  Piirlfcumnl,  1 970 

The  voluntary  informed  consent  of  the  subjects  used  In  this 
research  wn-.  i.otnined  in  iiccontamc  with  AFR  XO-JJ, 


useful  model  of  human  blood-pressure  response  to  (i 
stress  have  hitherto  not  been  available  An  efficient 
method  for  obtaining  a description  of  human  arterial 
oxygen  saturation  during  +G7.  stress  was  recently  re- 
ported (4):  ratios  of  Fourier  transforms  of  output 
oxygen  saturation  to  input  G stress  provided  frequency 
response  curves — i.c„  transfer  functions  describing  flu- 
general  mathematical  relationship  between  those  vari- 
ables. Despite  the  fact  that  the  biological  system  under 
lest  was  brashly  assumed  to  be  linear  and  stationary, 
surprisingly  consistent  descriptions  with  reasonable  pre- 
dictive ability  were  obtained.  By  employing  a method 
of  analysis  similar  to  that  reported,  we  sought  a work- 
able transfer  function  for  the  human  evc-lcvel  blood- 
pressure  response  to  +G/  stress. 

MATERIALS  AND  METHODS 

Three,  healthy,  male,  24,  25,  and  26  year-old  volun- 
teer members  of  the  USAFSAM  Acceleration  Stress 
Panel  provided  the  data  reported  herein.  The  left  radial 
artery  of  each  was  cannulatcd  with  a 5-em  Median* 

I X-gauge  catheter,  and  an  adequate  length  of  saline- 
filled,  noncompliant,  15-gaugc  Teflon*  tubing  connected 
the  arterial  catheter  to  a Stathnm  P37  pressure  trans- 
ducer in-  unfed  on  a Velcro*  headband  at  eye  level. 
A three-way  stopcock,  inserted  between  fhc  arterial 
catheter  and  the  tubing,  allowed  periodic  flushing  or  the 
catheter  with  heparinized  saline.  The  output  voltage  of 
the  pressure  transducer,  corresponding  to  instantaneous 
cyc-lcvel  blood  pressure,  was  recorded  on  magnetic  tape. 

Subjects  were  exposed  in  a seated  (13°  sent  hack 
tingle),  relaxed  condition— without  benefit  of  anti-G 
suit  or  straining  mnneuver— to  three  different  categories 
of  G stress  on  the  centrifuge:  I)  grnduol-onsct  runs 
(GORs),  in  which  the  G load  increased  linearly  with 
time  at  the  rate  of  0.067  G/s;  2)  rapid-onset  runs 
(RORs),  in  which  the  G load  rose  nt  about  L0  G/s 
to  a predetermined  ievel  and  remained  at  that  level  for 
cither  15  s or  1 min;  and  3)  simulated  aerial  combat 
maneuvers  (SACMs),  in  which  O loads  were  applied 
In  u varying  manner  over  a 100-s  period.  The  SACM  G 
profiles  were  of  eight  different  shapes,  each  profile 
having  been  produced  by  a random  function  generator 
with  bandwidth  set  to  mntch  the  capabilities  of  the 
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By  enscmhlc  averaging  transfer  functions  obtained 
from  single  input-output  pairs,  a mean  transfer  func- 
tion, 11(f),  derived  from  all  pertinent  data  in  a data 
set  (c.g,.  one  type  of  (i-stress  profile)  was  readily  ob- 
tained, Spectral  inconsistencies  (noise)  were  appreciably 
reduced  by  this  process.  Mean  responses  to  specific 
input  functions  were  obtained  by  averaging  the  transfer 
functions  generated  by  that  input  and  the  responses  to 
it,  and  transforming  to  the  time  domain  the  product 
of  the  average  transfer  function  and  the  Fourier  trans- 
form of  the  input;  i.c., 

6(f)  ^ ii(f)  Kf)  (Fq.  :> 

where  ()  (f)  is  the  Fourier  transform  of  the  mean  out- 
put, oft). 

Multiplying  a transfer  function  by  the  Fourier  trans- 
form of  a particular  input  of  interest  provided  the 
predicted  response  to  that  input,  based  on  the  given 
transfu  function,  In  oilier  words, 


UNAFSAs*  Human  C’entrfugc  (Fig.  I);  some  profiles 
were  augmented  or  clipped,  as  necessary,  to  challenge 
the  subject  visua'  functioning  without  rendering  them 
unconscious, 

Subjects  were  instructed  to  relax  completely  during 
a!!  lest  i ns.  and  to  release  a hand-held  enabling 
switch --(Ids  (ermnating  the  run — whenever  they  ex- 
perienced in*  visual  endpoint,  which  consisted  of  total 
loss  of  pe  pheral  vision  plus  significant  deterioration  of 
central  ve  :on.  Each  subject  was  exposed  to  the  follow- 
ing: one  < ;OR  to  visual  endpoint;  three  to  five  1 5-s 
RORs,  tin  final  one  eliciting  the  visual  endpoint;  at 
least  one  t )-s  ROD;  and  four  to  1 ! SACMs,  several  of 
which  elt.  ted  the  visual  endpoint  in  two  of  the  three 
subjects,  ' he  one  subject  who  sustained  only  four 
SACMs  t miniated  the  exocrimenl  early  because  of 
motion  si.  tncss. 

Disere!  , finite  Fourier  transforms  of  200-s  records 
of  input  i i stress)  and  output  (eye-level  arterial  blood 
pressure)  were  obtained  with  a Hewlett-Packard  545 IB 
Fourier  Analyzer,  !o  convert  the  instantaneous  blood- 
pressure  s ima!  to  mean  blood  pressure,  and  to  prevent 
spectral  a i.tsing  tv  the  transformed  data,  both  input 
and  output  signals  were  low-pass  filtered  with  a 0.5 -Hz, 
4X-dB/oenve  Butt  -rworth  filter.  The  paired  input  and 
output  waveforms  were  sampled  256  times  during  the 
200-s  aniuysis  period,  T;  i e„  At  » 0,78  s,  The  spectra 
obtained,  therefore,  had  a maximum  frequency  of 
0.04  Hz  At) 1 1 and  a frequency  resolution  of 

5,0  rnllz  (Af  » 1 /T).  A*  the  transient  time-domain 
signals  of  interest  began  ami  ended  well  within  the 
200-s  analysis  period,  the  octangular  time  window  used 
was  optimum  for  the  presen  application. 

The  transfer  function,  1(f),  between  a G-stress  Input, 
ift),  and  f lood-pressure  output,  o(t),  wits  obtained  by 
dividing  the  Fourier  transform  of  the  output,  0(f),  by 
that  of  the  input,  Iff);  i.c,, 
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O'(f)  - H ( f ) I'(f)  (Fq.  3) 

where  I '(f)  is  the  Fourier  transform  of  the  input  of 
interest,  i'(t),  and  O’(f)  is  the  Fourier  transform  of 
the  predicted  response,  o'(t). 

Because  the  power  spectra  of  the  SACMs  extended 
to  relatively  high  fiequcneies  without  tegular  spectrum 
nulls,  the  mean  transfer  function  obtained  by  enscmhlc 
averaging  the  transfer  functions  from  the  23  completed 
and  aborted  SACMs  was  examined  most  thoroughly  and 
with  the  greatest  expectations,  This  averaging  process 
gave  each  of  the  23  SACM  runs  equal  weight  in  the 
mean  transfer  function.  The  alternative  method  - 
averaging  each  subject's  SACM-gcncrated  transfer  func- 
tions and  then  averaging  the  three  subjects*  mean 
tmnsfer  functions — was  also  accomplished,  but  this 
process  theoretically  gave  undue  weight  to  the  noise 
associated  with  the  responses  of  the  subject  exposed  to 
only  four  SACMs,  I-ach  subject’s  individual  transfer 
functions,  as  well  ns  the  three  subjects'  mean  transfer 
functions,  were  sufficiently  similar  that  the  two  methods 
of  averaging  the  SACM-genernted  transfer  functions 
yielded  nearly  identical'  overall  mean  transfer  functions. 

Having  obtained  the  empirical  trnnsfer  function  re- 
lated eye-level  blood  pressure  to  +(iz.  stress  from  the 
ensemble  average  of  the  SACM-gcncrated  transfer  func- 
tions, we  elected  to  find  jevcrul  mathematically  well- 
behaved  substitutes  for  the  empirical  trnnsfer  function. 
The  analytic  transfer  functions  sought  were  of  three 
forms:  I ) single-zero,  double-pole;  2)  double-zero, 
double-pole;  and  3)  single-zero,  double-pole  with  a 
delayer,  Parameter  estimates  were  obtained  for  each 
analytic  form  by  minimizing  the  weighted  sum  of  the 
squared  error  of  the  fit  of  real  and  imaginary  terms  of 
the  analytic  form  to  the  first  40  frequencies  of  the 
empirical  trnnsfer  function,  ’She  weight  used  at  each 
frequency  was  inversely  proportional  to  the  variance  of 
the  data  comprising  the  empirical  function  at  that  fre- 
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Hg.  2.  Mean  Ittood-pressine 
transfer  functions  dcrivcil  front 
SACM  (i-slicss  runs  of  subjects 
Smith,  Mtntitlt.in,  ami  Carpenter. 
Him  50  frequencies  5 to  250  mil/ 
— displayed  at  5-mll/  Af,  Unedited 
data, 


qur  :y,  Marquardl’s  algorithm,  combining  the  Taylor 
series  expansion  ami  negative  gradient  methods,  was 
user  in  the  computations 

[•  tally,  predictions  of  eye-level  blood-pressure  rc- 
st'o’’  es  to  GOR,  ROR,  ami  SACM  (i  stresses,  based 
on  > te  empirical  and  various  analytic  transfer  functions, 
we:  compared  ills  the  mean  aelttal  responses  to  the 
san  (i-stress  inputs.  Such  comparisons  gave  visual 
ecu  me  of  the  tehtlive  merits  of  each  of  the  transfer 
furt'  sons  being  studied. 


R!  'IIS 


ach  subject's  mean  (i-to-blood-pressure  transfer 
fun  lion  obtained  from  SACM  Ci-strcss  runs  is 
shown  in  linear-scale,  polar  form,  in  Tig.  2.  Note  the 
similar  frequency  characteristics;  namely,  the  apparent 
resonance  peak  at  approximately  (>()  mHz;  the  tendency 
toward  — 1 H0°  phase  at  low  frequencies,  a phase  cross- 
over to  •*  I HO’  at  about  40  mHz,  and  a progressive 
dec-ease  in  phase  angle  at  higher  frequencies.  The  mean 
transfer  function  for  all  23  SACM  Ci-strcss  runs,  i.e, 

I 23 

(Iff)  « v ‘ H,(f) 

23  i ■**  I 

is  presented  in  Fig.  3.  The  frequency  charm,  .eristics  that 
were  seen  in  the  subjects’  mean  transfer  functions  arc 
again  evident,  as  is  the  improvement  in  the  signal-to- 
noisc  ratio  resulting  from  the  additional  averaging,  In 
this  overall  mean  transfer  function,  the  resonance  peak 
is  seen  to  be  at  fiO  mHz,  the  phase  crossover  occurs 
at  between  40  and  45  mHz,  and  a local  maximum  in 
the  pliasi  characteristic  is  noticeable  at  25  mHz.  In 
addition,  a 4-dIt/nctuvc  decline  in  magnitude  can  be 
estimated  for  frequencies  beyond  the  resonance.  Above 
200  mHz  no  useful  information  appears  in  the  frequency 
response  curve,  as  input  pmver  associated  with  the 
SACM  ('-force  profiles  was  virtually  nonexistent  above 


this  frequency.  The  first  40  terms  of  the  23-run  mean 
transfer  function  arc  given  in  Table  !. 

The  relatively  smooth  contours  of  the  above  empiric- 
ally determined  transfer  function  enticed  us  to  search 
for  analytic  functions  which  could  effectively  duplicate 
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crossover  occurs  between  50  and  55  mil/.  A con 
tirumif*  high-ficqucncy  decline  in  magnitude  is  not  pro- 
vided by  this  transfer  function,  although  a decline  com- 
paraidc  to  that  of  the  empirical  transfer  function  is 
evident  in  the  60-  to  200-mlf/  ratine.  After  the  phase 
crossover,  the  phase  angle  falls  and  rises  again;  Such 
behavior  is  not  necessarily  incongruous  with  that  of  the 
phase  characteristic  of  the  empirical  transfer  function 
since,  in  the  latter,  noise  may  be  obscuring  a slight 
upward  trend. 

Finally,  we  fitted  to  the  empirical  transfer  function 
an  analytic  function  consisting  of  one  zero,  two  poles, 
and  an  exponential  term: 


H(s)  - -IR. | 


1 + 7.66s 
I + 4,46s  + 7.63s1' 


» »<u« 

e tl'(|.  6) 


The  positive  exponent,  serving  to  augment  the  phase 
angle  in  proportion  to  frequency,  implies  anticipatory 
behavior  of  the  function:  our  observation  of  occasional 
rises  in  blood  pressure  prior  to  onset  of  Ci  stress 
prompted  us  to  evaluate  this  type  of  transfer  function. 
The  fitted  function  (Fig.  6)  displays  the  resonance 
peak,  phase  crossover,  and  other  features  of  the  em- 
pirical transfer  function  derived  from  the  23  SACM 
runs, 

We  noted  with  satisfaction  that  the  steady-state  char- 
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the  empirical  one.  First  to  be  evaluated  was’ a single- 
zero,  double-pole  transfer  function: 


H(s) 


! 6.2  — 


4-  1 2.0s 


! + 6.74s  + 7.11s* 


(F.q.  4) 


where  s is  the  conventional  complex  frequency  variable, 
j2rrf.  This  transfer  function,  presented  in  linear-scale, 
polar  form  in  F ig.  4,  shows  n resonance  peak  at  60 
mHz,  a local  maximum  in  phase  angle  at  15  mHz,  a 
phase  crossover  betwee  40  and  45  mHz,  and  high- 
frequency  magn  tude  tin  ',  phase-angle  reclines  similar 
to  those  of  the  empirical  transfer  function. 

'he  double-ze  o,  double-pole  transfer  f notion, 


H(s) 


18.3 


1 :•  7,1 8s  + 4.3{!  * 


I + 3.99s  + 7,9:  s* 


diq.  5) 


is  shown  in  Fig,  5,  Again,  resonance  is  seen  at  60 
mHz,  but  the  resonance  peak  is  sharp  tr  than  in  the 
si- -tie-zero,  double-pole  transfer  function.  The  local 
p1  ne-nnRle  maximum  is  seen  at  25  mHz  and  the  phase 
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FI*.  4,  Smite-zero,  double  xile  transfer  function  fitted 
empirical  transfer  function  In  Fig,  3.  Parameters  In  text. 
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acteristics  nf  t ic  second  two  fitted  analytic  transfer 
function'-,  came  out  to  be  - :K,7  and  — IK. I — not. too 

far  from  the  hcnrctical  . 2 mm  I Ip/CJ  discussed 
earlier, 

Hy  ‘akim*  fie  inverse  nttier  transform  of  the 
iransfe  funciio  for  ;t  systc  , one  obtains  its  impulse 
resporv  hit),  'he  iime-dom  o response  of  the  system 
to  a h reinp  fit  action  of  inf  ;!e  amplitude  and  infini- 
tesimal duratioi  The  inverv  ransform  of  the  complete 
empirii  at  traits!  r function  ■■  -s  too  corrupted  by  high- 
frequency  noisi  for  it  meaningful  impulse  response  to 
be  recognized,  and  invere  transforms  of  truncated 
version  , of  the  empirical  ttauslct  I unction  were  some- 
what distorted  cy  truncation  artifacts  (Gibbs’  phenome- 
non). ihe  invest  Fourier  iransfotms  of  the  analytic 
transfei  functions,  however,  appeared  as  reasonable  im- 
pulse responses,  fog,  7 is  the  mi;  ulsc  response  corrc- 
spondinn  to  the  double-zero  douKe-polc  iransfe r func- 
tion. t '•  invere  translorms  of  t!  ,■  other  tsvo  analytic 
Iransfe  functions  were  similar.  he  initial  negative- 
going  • lood-pr  ssure  response  to  the  theoretical  posi- 
tive ('  force  in  pulse  applied  at  ; — 0— is  readily  ap- 
parent as  is  the  zero-line  crossing  nt  4.0  s and  the  sub- 
seque-  reborn  ! peak  at  .0  s.  The  high-frceucncy 
nseill  turn  tit  the  impulse  response  is  ait  artifact  re- 
sultin' from  th  • unavo  lahle  truncation  of  the  infinite 
analy  it  transb  functio.-. 


FREQUENCY  imH/1 

f' if",  ft.  Single-zero,  double-pole,  espnnenli.it  transfei  function 
fitted  to  empirical  transfer  function. 


Since  Ihe  utility  of  a transfer  function  is  determined 
not  by  how  mnthemalicnlly-well  behaved  it  is  but  by  how 
well  it  predicts  responses,  we  compared  the  predicted 
responses  of  the  23-run  mean  empirical  and  the  three 
analytic  transfer  functions  to  the  subjects’  mean  actual 


0 5 10  14  20  25  30 

TIME  fiecondtl 


big,  7.  Inverse  Courier  transform  (impulse  response)  of 
double-zero,  double-pole  transfer  function  Alt  US  available 
frequencies— 3 in  M0  mil*— were  used  in  < (imputation. 
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I ii'.  I’reilicicil  ;iml  mean  aciuat  eyc-lcvcl  blood -pressure  re- 
spnnsc'  to  (!OR  (1  stress.  In  this  and  l her  following  figures,  top 
o.ue  i prediction  based  on  single-zero,  double-pole  transfer 
lime  no.,  seen  tul.  sinide-zeto,  double-pole,  exponential;  third, 
ifonhle-/eii).  double-pole;  loiiith,  2,1-run  mean  empirical.  I- if  tit 
ttace  is  mean  actual  i espouse  of  three  subjects -to  Ci  stress  in- 
ti,call'd  in  bottom  trace. 


l id.  Predicted  and  mean  actual  eye-level  blood-prcssuc  re- 
sponses to  aborted  ROR  (i  stress.  Sequence  as  m Rig.  8. 


responses  t<>  ,t  variety  of  G-stress  profiles.  The  predicted 
responses  to  a typical  CiOR  G-stress  profile  (Fig.  8), 
although  quite  similar  to  each  other  in  the  amount  of 
hlootl-pressure  drop  at  peak  G and  amount  of  over- 
shoot following  deceletation.  were  disappointingly  not 
very  well  matched  to  the  mean  actual  response  during 
die  period  of  expected  blood-pressure  rebound,  For- 
lunately,  the  reason  lor  the  discrepancy  lies  not  with 
the  transfer  functions  but  with  the  .subjects’  mean  re- 
sponse to  the  CiOR  G stress  only  one  response  to  this 
type  of  stress  was  obtained  from  each,  subject,  and  two 
of  the  three  responses  were  decidedly  atypical  when 
compared  with  blood-pressure  responses  to  GOT  G stress 
obtained  in  this  laboratory  on  oilier  occasions. 

More  sa  jsfaclory  was  the  comparison  of  the  various 
prediction  with  the  mean  actual  response  to  ROR  G 
'■tress.  Hot1  the  predicted  and  mean  actual  responses  to 
an  aborted  ROR  (Mg.  9),  bearing  understandable 
similarity  to  the  impulse  responses  of  the  analytic  trans- 
fer functions,  displayed  an  initial  blood-prcssurc  drop 
proportional  let  (I  In  d and  an  overshoot  occurring 
during  and  afier  deceleration,  The  prerun  rise  and  the 
postrun  oscillation  in  eresstue  seen  in  the  mean  actual 
response  were  not  pre  liclcd  by  cither  the  empirical  or 
the  analytic  transfer  functions,  A particularly  rigorous 
test  of  the  predictive  abilities  of  the  various  transfer 
functions  was  the  60-s  ROR  G profile  (Fig.  10),  since 
subtle  differences  in  the  responses  were  accentuated  by 
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Fig.  10.  Predicted  and  mean  actual  eye-level  blood-prenun 
response*  to  60-s  ROR  O stress,  Sequence  as  in  Fig,  a. 
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1 - j h II.  I'reihcteil  anil  mean  actual  eyc-'evc!  blood-pressure  re- 
sponses In  it  8 SAC  M (i  Mi  ass.  Scqu  .nice  as  in  Fig.  8.  To  avoid 
cnnlaminalmg  prediction  woh  actual  response  data,  data  from 
caclt  subject's  itH  SACM  iun  were  subtracted  from  empirical 
transfer  ftuunon  pnor  to  making  prediction. 


the  long,  straight  lines  of  the  forcing  function.  The 
mean  actual  response  to  ROR  G stress  showed  a numbe; 
of  characteristics  coinciding  with  expectations  based  on 
prior  experience:  1)  an  anticipatory  rise  in  blood  pres- 
sure prior  to  onset  of  G stress;  2}  a rapid  fall  in  pressure 
at  onset  of  G stress;  3)  a partial  recovery  of  pressure 
several  seconds  after  the  G load  had  stabilized;  4) 
a secondary  gradual  decline  in  pressure  following  the 
peak  of  the  partial  recovery;  5)  a rapid  rise  in  pressure 
at  the  offset  of  G stress;  6)  an  overshoot  in  pressure 
coincident  with  termination  of  the  G load;  and  7)  a 
small  unde i shoot  immediately  following  the  terminal 
overshoot.  Whereas  the  empirical  and  analytic  transfer 
functions  all  predicted  the  pressure  drop  at  onset,  pres- 
sure rise  at  offset,  partial  recovery  at  stabilization,  and 
overshoot  at  termination  of  G stress,  they  met  'with 
varying  degrees  of  suet  ess  in  predicting  other  com- 
ponents of  the  waveform.  Not  one  recreated  the  antici- 
pattny  rise  in  eye-level  blood  pressure.  Only  the  em- 
pirical transfer  function  predicted  the  gradual  decline  in 
pressure  during  the  constant  G stress,  and  only  the 
empirical  provided  a noticeable  undershoot  following 
the  terminal  overshoot.  It  is  worthy  of  note  that  the 
predie  id  responses  to  ROR  G stress  had  the  same 
gencr  characteristics  a the  classic  figures  picsented  by 
Wood  and  others  (11.16);  in  fact,  the  greatest  con- 
gruity  vith  Wood's  data  was  obtained  with  the  empirical 
transf*.  r function,  in  that  it  effectively  reproduced  the 


G(») 


R(s> 


PM 


Fig.  !2.  Block  diagram  of  eye-level  blood-prewure  control 
system  for  relaxed  humans  exposed  to  +(iz  stress.  A(»),  B(s). 
t'(s),  and  D(s)  are  the  hydraulic,  neurocirculatory  effector, 
baroceptor,  nnd  vestibulocerebellar  components  of  the  control 
system,  respectively,  R(«)  ix  the  reference  signal,  l*(»)  it  eye- 
level  blood  pressure,  and  G(x)  is  G-streis  perturbation. 


gradual  decline  in  pressure  during  the  latter  portion  of 
the  plateau.  Of  the  analytic  transfer  functions,  both  the 
double-zero  and  tllc_cxponential  forms  yielded  more  ac- 
curate predictions  than  did  the  single-zero  form,  which 
provided  too  smooth  a transition  into  the  partial  re- 
covery phase  during  G stress  and  loo  gradual  a poslrun 
return  to  mean  pressure. 

By  averaging  the  individual  transfer  functions  ob- 
tained from  all  SACM  G-stress  runs,  except  that  SACM 
run  for  which  a prediction  is  to  be  made,  a mean 
transfer  function  uninfluenced  by  the  test  SACM  can 
be  obtained.  Such  transfer  functions  were  used  to  make 
predictions  of  eye-level  blood-ptessurc  response  to  the 
several  different  forms  of  SACM  G stress.  Fig.  1 1 
shows  the  mean  actual  response  to  the  #8  SACM,  as 
well  as  responses  to  this  SACM  predicted  by  the  three 
analytic  transfer  functions  and  by  the  empirical  transfer 
function  with  #8  SACM  data  excluded.  Although  it  was 
difficult  to  evaluate  visually  the  relative  accuracy  of 
predicted  responses  to  so  complex  a waveform  as  the 
#8  SACM  profile,  such  characteristics  as  the  relative 
depths  of  the  first  two  negative  peaks,  the  relative 
heights  of  the  highest  two  positive  peaks,  and  the  shape 
of  tlte  recovery  from  the  terminal  rebound  enabled  us 
to  observe  that  the  empirical  transfer  function  was  the 
most  accurate, 

To  determine  the  goodness  of  fit  of  the  three  analytic 
transfer  functions  to  the  23-run  mean  empirical  transfer 
function,  the  standard  error  of  estimate  was  obtained 
for  each.  This  figure  of  merit  was  4.02  for  the  single- 
zero,  double-pole  mode!;  3,34  for  the  single-aero, 
double-pole,  exponential  model;  and  3.26  for  the 
double-zero,  double-pole  model,  The  doublc-icro, 
double-pole  transfer  function  thus  emerged  as  the  best 
of  the  analytic  functions  tested  in  that  it  most  closely 
matched  the  empirical  transfer  function. 
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discussion 

We  recognize  t ho  1 1 in i l:il ions  of  us j nj*  linear  systems 
analysis  techniques  on  notoriously  nonlinear,  nonsla- 
tionaty,  biological  systems.  The  Mood-pressure  control 
system  was  rensonaMy  traetahle  to  Fourier  tinttlysis, 
however,  since  well-behaved  frequency  response  func- 
tions could  he  recognized  ttttd  substantial  predictive 
ability  of  derived  transfer  functions  could  he  demon- 
Minted,  With  regard  to  selection  ol  the  most  appropriate 
transfer  function,  it  appears  that  the  fit st  40  terms  of 
the  empirical  transfer  function  would  he  the  logical 
clause  if  one  requires  greatest  predietive  accuracy  and 
closest  ties  to  experimental  data.  If  concise  mathematical 
expression  is  of  greater  importance,  the  double-zero, 
double- pole  ttansier  Itmction  should  he  sufficient.  Use 
ol  the  single  zeto,  dnuhle-polc,  exponential  transfer 
hmetion  would  not  he  justilted  since  no  advantage  of 
the  exponential  term  over  an  additional  zero  is  demon- 
strable. 

A simplified  block  diagram  for  the  eye-level  blood- 
pressure  control  system  in  relaxed  lutmans  exposed  to 
U stress  is  shown  in  Fie.  12.  The  (i-to-pressure  transfer 
function,  I Its),  is  thus  described  in  litis  scheme  as 


„ A ( s ) + M(s)Ws) 

Ilfs)  

I + H(  s)C'(s) 


(F-q.  7) 


where  Afs)  repiesents  the  direct  effects  of  the  hydro- 
static load  on  the  Mm>d  and  body,  H(s)  includes  the 
various  neural,  hormonal,  and  mechanical  effectors  that 
translate  control  signals  into  blood  pressure  (primarily 
via  changes  in  heat:  rate,  venous  compliance,  and 
arterial  resistance).  and  ('(s)  represents  the  dynamics 
of  the  several  baroceptor  feedback  mechanisms.  The 
effects  (>f  the  vestibulocerebellar  inputs  to  the  cardio- 
vascular control  system  are  included  in  D(s).  "1  he  hy- 
drosltttie  column  ellcc;  per  xc  is  frequency-independent, 
so  th  t at  very  low  bequencies  of  G stress,  the  hydro- 
static -011111111  effect  predominates;  at  higher  frequencies, 
the  ! equeney-dependent  properties  of  compliant  body 
st  i tic.  j res,  vestibulocerebellar  pathways,  cardiox  asetdar 
eflec  >r  mechanisms,  ttttd  the  baroceptor  feedback  loops 
beco:  c increasingly  important.  The  properties  < f pure- 
ly sivc  elements,  i.e.,  the  elasticity  of  the  vascular 
tree  tul  the  inertia  of  the  enc'osed  blood,  provide  a 
potential  for  resonance,  which  is  the  reason  Afs)  is 
denoted  as  a complex  variable  in  the  ( i-to-bloc id-pres- 
sure ransfer  function.  The  vestibulocerebellar  eftegts  on 
bloo  ! pressine  are  likely  proportional-plus-deriv;  tive  in 
font  , since  the  vestibulocerebellar  mechanisms  appear 
to  participate  mainly  in  initiating,  but  also  in  sustaining, 
the  orthostatic  reflexes  (h.  This  implies  that  Dfs)  in 
the  i-to-blood-pressurc  transfer  function  is  u term  no 
simpler  than  K,(l  •)-  K;  •).  The  open-loop  transfer 
chat  acteristies  of  the  carotid  sinus  baroceptor  system 
have  been  studied  in  animals  (K,I0,I3),  and  the  evi- 
dence indicates,  at  the  very  least,  an  integrating  function 
of  ti  e form  K.,/(  I -t-  Kts)  where  the  time  constant  Kt  is 
on  t re  order  of  20  or  30  s,  The  linear  portion  of  one 
reported  model  (10)  based  on  the  input-output  be- 


havior of  the  canine  carotid  sinus  reflex  is  of  the  form 
K;,/(l  T K,s  + K.,s-):  interestingly,  this  double-pole 
model  has  a natural  frequency  of  42  mil/,  which  is 
not  far  from  the  resonance  point  of  the  empirical  trans- 
fer function  obtained  in  the  present  experiment.  Another 
such  model  (13)  even  suggests  a single  zero,  double- 
pole  loop  transfer  characteristic.  From  such  evidence 
we  should  certainly  expect  that  the  combined  transfer 
function,  incorporating  the  elastic  elements,  the  vesti 
htiloeercbcllar  pathway,  and  the  baroceptor  feedback 
loop,  contains  at  least  one  zero  and  two  or  more  poles. 

I he  evidence  for  additional  zeros  ami  poles  at  fre- 
quencies above  1.0  11/  (X)  must  he  taken  into  account 
in  the  inevitable  complete  model  of  the  blond -pressure 
control  system,  but  is  not  particularly  relevant  to  our 
discussion  of  blood-pressure  responses  to  low-frequency 
G stress. 

Knapp  ci  ul.  (7)  recently  plotted  the  frequency  ie- 
spouse  of  the  canine  cardiovascular  system  as  a whole. 
Using  ±2.0  Ci/  sinusoidal  input  stress  tit  discrete  fre- 
quencies from  0.001  to  1.5  II/..  they  obtained  the  G- 
lo-blood-pressurc  transfer  characteristic;  maximum  gain 
was  observed  at  30  to  60  ini  Jz,  in  rough  agreement 
with  on r own  data_on  the  human  response.  They  also 
reported  that  the  frequency  response  of  dogs  with  heart 
rate  kept  constant  and  peripheral  vascular  response 
blocked  pharmacologically  displayed  a maximum  be- 
tween 40  and  60  mHz,  confirming  the  complex  form  for 
the  purely  hydraulic  segment  of  the  blood-pressure  con- 
trol system.  Also  of  interest  were  their  observations 
that  the  influence  of  peripheral  mechanisms  was  present 
up  It)  30  mHz,  while  heart-rate  responses  were  most  ap- 
parent in  the  20-  to  100-mHz  range,  thus  providing 
additional  evidence  that  the  net  G to-blood-prcssute 
transfer  function  consists  of  multiple  poles  am!  zeros 

Such  data  as  are  presented  herein  for  relaxed  humans 
and  by  Knapp  ci  at,  for  dogs  arc  the  type  prerequisite 
for  the  development  of  models  like  that  of  Green  and 
Miller  (5)  and  that  of  Witte  cl  al.  (15).  We  hope  that 
those  who  develop  models  of  cardiovascular  function 
during  G stress  will  decompose  the  overall  transfer 
functions  presented  to  obtain  estimates  of  the  transfer 
characteristics  of  the  various  segments  of  the  Mood- 
pressure  control  system.  We  also  hope  these  data  will 
be  useful  to  modelers  of  pilot  performance  during  aerial 
combat  maneuvering.  Required  for  such  use.  however, 
is  a proper  description  of  the  changes  in  cardiovascular 
dynamics  resulting  from  straining  maneuvers  and  the 
use  of  unti-G  suits,  valves,  and  seats.  Since  struining 
maneuvers  and  the  operation  of  most  anti-G  protective 
equipment  arc  highly  nonlinear,  it  would  be  neccssatx 
to  model  these  effects  separately,  then  add  them  to  the 
basic  blood-pressure  control  system,  rather  than  to  try  to 
model  as  a linear  system  the  overall  mechanism  of  the 
protected  cardiovascular  response  to  G stress. 

CONCLUSION 

Dala  describing  the  cye-lcvcl  blood-pressure  rcs|H>nsc 
to  4-0*  stress  in  relaxed  humans  have  been  presented  as 
an  empirical  G-to-blood-prcssure  transfer  function.  This 
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transfer  fun.  lion  covers  the  frequency  spectrum  from 
•<i  ,u  ~(>0:  diM/.,  and  exhibits  substantial  predictive 
ability  in  Mu  t .()-•!.>  (I,  0-200  s ranges.  Three  simple 
m.ithemalica  models  of  the  (i-stressed  blood-pressure 
eontinl  system  were  also  tested  for  predictive  ability, 
and  were  found  moderately  valid.  Of  these,  the  double- 
zero,  double  pole  transfer  function  most  closely  fits  the 
empirical  transfer  function,  and  appears  to  predict  most 
accurately , 
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